Charge-density waves in one-dimensional Hubbard superlattices 
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We study the formation of charge density waves (CDW's) in one- dimensional Hubbard superlattices, 
modeled by a repeated pattern of repulsive (U > 0) and free (U = 0) sites. By means of Lanczos 
diagonalizations for the ground state, we calculate the charge structure factor. Our results show that 
while the superlattice structure affects the modulation of the charge density waves, the periodicity 
can still be predicted through an effective density. We also show that, for a fixed repulsive layer 
thickness, the periodicity of the CDW is an oscillatory function of the free layer thickness. 

PACS: 71.45.Lr, 73.21.Cd, 7f.27.+a, 71.10.-w, 73.20.Mf, 72.15.Nj. 



Charge-density waves (CDW's) are present in a va- 
riety of homogeneous strongly correlated electron sys- 
tems, especially some quasi-one-dimensional organic 
conductors .Orel An interesting property of CDW systems 
is their non-ohmic behavior: CDW's incommensurate 
with the underlying lattice are pinned by impurities, so 
that only voltages above a certain threshold are able to 
produce observable currents. Just above this critical volt- 
age, the current rises-sharply, indicating a huge increase 
in the conductivity.Effi 

Thus, one of the key issues in the study of CDW's 
is to determine the period of charge modulation in the 
ground state. From the theoretical point of view, this 
has been addressed with the aid of various models, but 
care must be taken even in the simplest case of the one- 
dimensional Hubbard model. Indeed, according to the 
Luttinger liquid (LL) description|j1i3 the large-distance 
behavior of the charge density correlation function in the 
ground state is given by 
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where the amplitudes A\ and A2, and the exponent K p 
are interaction- and density-dependent parameters; f< 
the Hubbard model with repulsive interactions one hasl 
f < K p < 1, so that charge correlations are expected to 
be dominated by the 2kp term, where 2kp — np is the 
Fermi wave vector for a density p < 1 of free electrons 
on a periodic lattice; if p > 1, then 2kp = ?r(2 — p). By 
contrast, numerical analyses predict that CDW correla- 
tions are dominated by .the AkF term, for strong enough 
on-site repulsion [/.EirEj These two scenarios-ean only be 
reconciled if A\ vanishes above some t/*(p).l!3 

On the other hand, the development of diverse and 
very accurate deposition techniques over the last two 
decades has generated a whole new class of materials, 
generically called heterostructures. Systems made up of 
very thin layers - in some cases a few atoms thick - of ma- 
terials with different properties have been grown; if the 
layered structure is periodic, it is referred to as a super- 



lattice (SL). Interesting new phenomena have emerged 
from these heterostructures, demanding novel theoreti- 
cal mechanisms to explain the observed data, and giving 
rise to improved devices for technological applications. 

We have recently considered a model for strongly parre- 
lated electrons in one-dimensional superlattices,E3ii3 de- 
fined by the Hamiltonian 
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where, in standard notation, i runs over the sites of a 
one-dimensional lattice, c ia (ci a ) creates (annihilates) a 
fcrmion at site i in the spin state a =f or |, and 
n>i = + riii, with rii a = c ia .Ci a ; the on-site Coulomb 
repulsion is taken to be site-dependent: Ui = U > 0, for 
sites within the repulsive layers, and Ui = otherwise. 
The repulsive and free layers have Ljj and Lq sites, re- 
spectively; the 'aspect ratio' is defined as i = Ljj/Lq. 
The extension of this model to a Luttingerdiauid super- 
lattice (LLSL) has been discussed recently.ESEZI 

The magnetic properties of model (||) turned out to be 
quite different r Jxorn those of the corresponding homoge- 
neous system;E3llj also, the electronic density at which 
a charge gap opens in these, S i Lj s can be fine-tuned by 
an appropriate choice of £.L3'E3llZl In view of this, the be- 
havior of CDW's on a superlattice structure should also 
be of interest. In particular, a crucial question is how 
to identify and parametrize the dominant charge-density 
correlations in the case of these Hubbard superlattices 
(HSL's). Our purpose here is to address these points. 

We use the Lanczos algorithmic c!J to determine the 
ground state of Eq. (||) , for finite lattices of N s sites with 
periodic boundary conditions, in the subspace of fixed 
particle-density (canonical ensemble) p — N e /N s , where 
N e is the total number of electrons. The signature of a 
CDW instability is a cusp at q = q* , in the charge-density 
structure factor, 
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FIG. 1. Charge structure factor for p — 11/6 and N a = 24, 
for a homogeneous chain (a), and for a superlattice with 
Lu — Lo — 1 (b). The curves are labelled by the values of 
the on-site coupling: U = (open circles), U — 2 (squares), 
(7 = 4 (circles), U = 12 (down triangles), and U = 20 (dia- 
monds). Successive vertical shifts of 0.05 have been imposed 
on the curves for clarity. 



with (QiQj) = {ipalniUjlipa} - (ipo\ni\ipa){^a\nj\ipo), 
where \ipo) is the ground state and N c is the number of 
periodic cells, N c — N s /Nb, for a basis with Nf, = Ljj+Lq 
sites. 

In order to study the behavior of C(q) in the widest 
possible range of values of Lq, we focus our discussion 
on the case of 24-site lattices, with electron density p — 
11/6. Figure [l] compares the charge structure factor for 
the homogeneous chain with that for a SL with Lu — 
1 and Lq = 1. From Fig. |l|(a) one sees that the 2kp 
cusp in the free (U = Ql homogeneous case is displaced 
to as U increases.E3 The SL structure changes the 
periodicity of the CDW as shown in Fig. [l](b): the cusp 
is displaced relative to the position for the homogeneous 
chain. Similarly to the homogeneous system, the cusps 
sharpen as U increases; their positions, however, remain 
locked at 27r/3 for that SL configuration and electronic 
density. 

Before discussing the location of the cusps, let us check 
on the role of finite-size effects. Recall that the appropri- 
ate finite-size parameter is the number of repeating units, 
N c . Figure H shows C(q) for the above mentioned SL and 
for two different lattice sizes, namely N s — 12 (N c = 6) 
and N s = 24 (JV C = 12). One sees that for U = 4 the 
cusp is sharpened as one goes from N c — 6 to N c = 12; 
the cusp positions, however, do not change as the size 
increases. For large couplings (U = 30 in the figure), the 
curves lie on top of each other, so that finite-size round- 
ings are not noticeable. Thus we can rule out finite size 
effects as playing any crucial role in determining the cusp 
positions £3 
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FIG. 2. Charge structure factor for a superlattice with 
Lu = Lo = 1, p = 11/6, and both (7 = 4 (circles) and 
(7 = 30 (squares). Data corresponding to N s = 24 (N s = 12) 
are denoted by filled (open) symbols, joined by full (dotted) 
lines. For clarity, successive vertical shifts of 0.05 have been 
imposed on the curves corresponding to different values of (7. 



The main feature determining the periodicity of the 
CDW can be identified through a systematic study of 
q* as a function of layer thickness. Figure || shows the 
charge structure factor for a fixed length of the repulsive 
layer (Ljj = 1), and for different lengths of the free layer: 
Lo = 1, 2, 3, and 5. We see that the CDW modulation 
depends on the superlattice: q* fix = 2/3, 1, and 2/3 
for Lq = 1, 2, and 3, respectively; the vanishing C(q) 
obtained for Lq = 5 is associated with q* = 0. 
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FIG. 3. Charge structure factor for p = 11/6, (7 = 12, 
7V S = 24, Ljj — 1, and for different free layer lengths: (a) 
Lo = 1; (b) Lo = 2; (c) L = 3; (d) L = 5. 



These cusp positions can be analysed in a strong cou- 
pling scenario. We first recall that there is a special den- 
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sity, pi , corresponding to each free site being doubly oc- 
cupied, while.-each repulsive site is singly occupied; it can 
be written aalil 

At this density a MotUHubbard gap opens and the sys- 
tem is an insulator;E-M'LD that is, the transport proper- 
ties of a SL at pi should be-similar to that for a homo- 
geneous chain at-half-fillingjia which includes the break- 
down of CDW'sH This latter point is nicely illustrated 
by the vanishing of C (q) for all q for the SL with Ljj = 1 
and Lq — 5, as shown in Fig. ||(d): pj is exactly 11/6 for 
that SL. For densities above pi, all free sites remain dou- 
bly occupied, while_lhe amount of charge on the repulsive 
sites can fluctuate]!!] Charge correlations between differ- 
ent cells are therefore dominated by the way in which 
electrons are distributed within each repulsive layer. Ac- 
cordingly, we define the effective cell density (i.e., the 
number of active electrons per unit cell) , 

Pes = p(L + Lu) - 2L , (5) 

where p is the overall density. If one similarly defines an 
effective Fermi wave vector, k F , from 

2fc> = ^(2-p cff ), (6) 

then the cusp in C{q) is located simply at q* = 4k F , in 
a way completely analogous to that of the homogeneous 
system. In Fig. |l](b), 4k* F thus calculated is indicated by 
an arrow, and it is clear that it coincides with the cusp 
positions. By the same token, for all SL's shown in Fig. 
[|, the position of the cusps also perfectly match the 4k F 
as obtained from Eq. (|J). 




FIG. 4. q* as a function of Lq for iVs=24, U = 12, 
p — 11/6, and for different repulsive layer lenghts: (a) Lu = 1, 
and (b) Lu = 2. 



It is now illustrative to discuss how the periodicity of 
the CDW depends on the width of the free layer as the 
repulsive layer width is kept fixed, as done in Fig. ^ where 
q* is shown as a function of Lq for both Lu = 1 and Ljj — 
2. First, we see that the strong coupling prediction, q* = 
4k p, with kp given by Eq. (|6|), is satisfied for any U > 
0: it agrees perfectly well with the data points (circles) 
obtained from Lanczos diagonalization of different SL's 
with U — 12. And, second, we see that q* is periodic 
with Lq, for a fixed Ljj. The period of oscillation ALo 
can be determined by imposing q*(Lo) = q*{L$ + ALo), 
which yields 

ALo = (2-p)-* = JL (7) 

where p is the overall density. 

A similar oscillatory behavior has hfien found for the 
periodicity of SDW's on Hubbard SL's,E£l which was iden- 
tified with the exchange oscillation in magnetic multilay- 
ers; this latter feature is, in turn, intimately connected 
with the Giant Magnetoresistance (GMR) effect (see e.g., 
Ref. |2^ and references therein). Thus, the above results 
indicate that one can also fine-tune a given charge dis- 
tribution by a suitable choice of the length of the spacer 
(i.e., free) material. 

As a final remark, we should mention that experimen- 
tal realizations of (variations of) the model discussed 
here have recently appeared in the literature. Indeed, 
by growing La layers of semimetallic (narrow gap semi- 
conductor) TiSe2 in succession to Lb layers of metal- 
lic/superconducting NbSe2 - beth compounds undergo 
CDW transitions - Noh et alEn were able to grow 'su- 
perfilms' 1000 A thick, with La and Lb in the range 1 
to 24 layers. The superconducting critical temperature 
and the temperature dependence of the conductivity were 
then obtained as functions of the SL structure. 

In summary, we have considered the formation of 
CDW's in one-dimensional Hubbard superlattices. We 
have established that the SL structure affects the charge 
distribution in relation to the homogeneous lattices, 
though one can still predict the new periodicity through 
an effective cell density. It has also been established that 
the periodicity of the CDW is an oscillatory function of 
the spacer thickness, for a fixed repulsive layer thickness. 
As a consequence, for a given (overall) electronic den- 
sity, one can fine-tune a desired charge distribution by 
a suitable choice of the length of the spacer (i.e., free) 
material. We hope the present work stimulates further 
investigations on the superlattices made up of transition 
metal dichalcogenides, focusing on the CDW periodicity 
along the direction of growth. 
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